Non-Universal Spectra of Ultra-High Energy Cosmic Ray Primaries and Secondaries 
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Analytical calculations of extra-galactic cosmic ray spectra above ~ 10 17 eV are often performed 
assuming continuous source distributions, giving rise to spectra that depend little on the propagation 
mode, be it rectilinear or diffusive. We perform trajectory simulations for proton primaries in the 
probably more realistic case of discrete sources with a density of ~ 10 -5 Mpc -3 . We find two 
considerable non-universal effects that depend on source distributions and magnetic fields: First, 
the primary extra-galactic cosmic ray flux can become strongly suppressed below a few f 18 eV due 
to partial confinement in magnetic fields surrounding sources. Second, the secondary photon to 
primary cosmic ray flux ratio between ~ 3 x f 18 eV and ~ 10 20 eV decreases with decreasing source 
density and increasing magnetization. As a consequence, in acceleration scenarios for the origin of 
highest energy cosmic rays the fraction of secondary photons may be difficult to detect even for 
experiments such as Pierre Auger. The cosmogenic neutrino flux does not significantly depend on 
source density and magnetization. 
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I. INTRODUCTION 

A major unresolved aspect of ultra-high energy cos- 
mic ray (UHECR) physics [l[ is their composition and 
above which energy the flux is dominated by extragalac- 
tic sources. Above ~ 10 17 eV the chemical composition 
is basically unknown ||. Around 10 18 eV the situation is 
particularly inconclusive as HiRes [|| and HiRes-MIA Q 
data suggest a light (proton dominated) composition, 
whereas other experiments indicate a heavy composi- 
tion @. 

As a consequence, there are currently two different sce- 
narios: The "standard" one, where a transition from a 
steeper, galactic heavy component to a flatter, extra- 
galactic component dominated by protons takes place at 
the ankle at ~ 5 x 10 18 eV, see, e.g., Ref. @, and a more 
recent one suggesting that this transition actually takes 
place at lower energies, namely around the "second knee" 
at ~ 4 x 10 17 eV. 

This second scenario in which extragalactic protons 
dominate down to the second knee has the following con- 
sequences: First, since the observed spectrum above the 
second knee is quite steep, oc E~ 3 3 , the extragalactic 
proton flux has to cut off below ~ 4 x 10 17 eV. This can be 
explained as a magnetic horizon effect: Protons from cos- 
mological distances cannot reach the observer any more 
within a Hubble time due to diffusion in large scale mag- 
netic fields 0, HUES [UEl. Second ' the ankle would 
have to be interpreted as du e to pair production of the 
extragalactic protons 



pointed out recently 




. In particular, it has been 
that this model cannot af- 
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ford injection of a significant heavy component above the 
ankle whose secondary photo-disintegration products in 
the form of intermediate mass nuclei would produce a 
bump around the ankle. Injection of a mixed composition 
would also require a harder injection spectrum oc E~ a 
with a ~ 2.2, as opposed to the extragalactic ankle sce- 
nario with pure protons [14j which requires an injection 
spectrum oc E~ 2 6 . 

In the proton dominated low-cross over scenario the 
resulting spectra are surprisingly insensitive to details 
such as extra-galactic magnetic fields (EGMF) and ac- 
tual source distribution, as long as the distances between 
UHECR sources are much smaller than the energy loss 
and diffusion lengths [HI [l7| • The spectra are then equal 
to the spectra resulting from a homogeneous source dis- 
tribution which is therefore called "universal" . Continu- 
ously distributed sources fulfill this condition. 

Actual UHECR sources are, however, likely dis- 
crete with a relatively small density of order ~ 
10 ~ 5 Mpc -3 These values are motivated by the den- 
sity of candidate sources such as active galaxies and also 
by hints of UHECR clusters. In addition, there likely 
is a strongly structured EGMF. In the present paper we 
study how a strongly structured Universe can cause de- 
viations from the universal spectra in the case where ex- 
tragalactic protons dominate the observed flux down to 
- 10 18 eV. 

We also include in this study secondary neutrinos and 
7— rays, produced by interactions of the primary nucle- 
ons with the cosmic microwave background (CMB) and 
the infrared (IR) background. The relevant reactions 
which also govern energy loss of the primary nucleons are 
pion production above the "GZK threshold" [l9| and, for 
7— rays, by pair production of protons. There is a strong 
motivation to study predictions for the fraction of 7— rays 
in the UHECR flux from recent experimental upper lim- 
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its @, HI, H US El that will further improve in the 
near future, especially from the Pierre Auger project [25j . 
Even secondary neutrino fluxes may be detectable in the 
not too far future [2o| . 

In section 2 we describe our simulations, in section 
3 we discuss results for various degrees of structure in 
the source distribution and magnetic field scenarios. In 
section 4 we discuss uncertainties in the 7— ray to charged 
primary flux ratio predicted in bottom-up scenarios and 
we conclude in section 4. We use natural units, Ti = c = 
k = 1, throughout the paper. 



duced, are only subject to redshift. 

The concordance cosmology is used for which, as- 
suming a flat Universe, the Hubble rate H(z) at red- 
shift z in the matter dominated regime, z < 10 3 , is 

given by H(z) = H [O m (l + z) 3 + fi A ] 1/2 ■ We use 
the standard values being Q m — 0.3, Cl\ — 0.7, and 
H = h 100 km s" 1 Mpc" 1 with h = 0.72. 

III. RESULTS 



II. NUMERICAL SIMULATIONS 

For the trajectorysimulations we use the public code 
package CRPropa [27J, |28[. It follows nucleon trajecto- 
ries in 3 dimensions in arbitrarily structured environ- 
ments and takes into account secondary neutrinos and 
electromagnetic (EM) cascades produced by nucleon in- 
teractions with the CMB and IR backgrounds. Note that 
trajectory simulations allow to treat both the rectilin- 
ear and diffusive regime without approximation and are 
thus more general than using the diffusion approxima- 
tion as adopted in Ref. [HI, E3] which in addition was re- 
stricted to a homogeneous EGMF. We will also compare 
the general case to the isotropic situation in the absence 
of UHECR deflection and source structure which can be 
simulated with the one-dimensional version of CRPropa 
in which primary and secondary particles just propagate 
along strai ght lines. This case has been considered, e.g., 
in Ref. [IS|,l30|. 

We inject protons with a spectrum E~ 7 up to 10 21 eV 
and adjust the injection index 7 to the data. We also 
consider scenarios where the injection power per comov- 
ing volume evolves as a power law oc (1 + z) m , where m 
is often called the "bright phase index" . We integrate up 
to redshift z = 3. 

For nucleons we take into account photo-pion produc- 
tion, and pair production on the CMB and IR back- 
ground, as well as redshift, and deflection in the cases 
where we consider an EGMF. 

The 7— ray interaction length above ~ 10 19 eV strongly 
depends on the density of the universal radio back- 
ground (URB) [3l|. We use the minimal estimate for 
the URB [32[ ■ For the IR background we use the model 
of Primack et al. [33| . The IR background is mostly rel- 
evant for pion production for relatively steep proton in- 
jection spectra, and for photon attenuation in the TeV 
regime. Since we are mostly interested in energies above 
~ 10 17 eV, for our purposes the IR background is less rel- 
evant for EM cascade propagation. Between ~ 10 17 eV 
and ~ 10 21 eV, the 7— ray attenuation length is smaller 
than 10 Mpc, and thus redshift effects are negligible 
at these energies. Furthermore, since we are interested 
in diffuse fluxes, we can treat the EM cascades as 1- 
dimensional. 

Neutrinos propagate along straight lines and, once pro- 
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FIG. 1: Upper panel: Fluxes of primary nucleons (black), 
secondary photons (red) and neutrinos per flavor (blue) for 
a homogeneous source distribution injecting an E~ 2 ' 6 proton 
spectrum up to 10 21 eV and up to redshift z = 3. The dashed 
lines are in absence of redshift evolution, the solid lines are 
for a bright phase index m = 3. No magnetic fields are taken 
into account and fluxes have been computed with the one 
dimensional version of CRPropa. AGASA data [34| are shown 
as triangels, HiRes data [35[ as crosses. Lower panel: Same 
over a wider energy range. Also shown is the diffuse 7— ray 
flux that has been observed by EGRET 

In Fig. [1] we show nucleon, 7— ray and neutrino fluxes 
for a Universe filled with homogeneously distributed 
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sources, with and without source evolution. As pointed 
out in Ref . [HI, [ItJ , the HiRes and AGASA spectra can 
be made to give a consistent ankle structure, coinci- 
dent with the theoretically predicted ankle position at 
~ 5 x 10 18 eV, by multiplying the HiRes energies with 
a factor 1.2 and the AGASA energies with a factor 0.9. 
This is best seen by multiplying the differential spectrum 
with E 3 , as shown in Fig. O The best fit to the observed 
flux down to ~ 2 x 10 17 eV then gives an injection spectral 
index of 2.6 < 7 < 2.7, where the lower value corresponds 
to to ~ 3 and the upper value to to ~ 0. 
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FIG. 2: Same as Fig. [Tj upper panel, but with the energies 
measured by HiRes and AGASA multiplied by a factor of 1.2 
and 0.9, respectively, and showing the differential nucleon flux 
multiplied by E s . 

As the lower panel of Fig. Q]of shows, the 7— ray flux at 
GeV energies can reach about ten per cent of the diffuse 
background observed by EGRET [36| in this scenario, 
particularly for significant source evolution. The energy 
in GeV 7— rays is comparable to the energy in primary 
cosmic rays at ~ 2 x 10 18 eV, where the proton energy 
attenuation length due to pair production becomes com- 
parable to the Hubble radius [31(. A diffuse 7— ray flux 
at this level should be testable in the near future [37j |. 

As in our earlier work in Refs. [38| . for the scenar- 
ios with structured source distributions and/or magnetic 
fields we will in the following use the unconstrained large 
scale structure (LSS) simulations based on Refs. [39l.l4oT|. 
A cross-section through the baryon density and mag- 
netic field strength of the LSS simulation box is shown 
in Fig. [3] The EGMF in these simulations have been 
evolved passively and normalized to ~ /iGauss in the 
centers of galaxy clusters. These EGMF are highly 
structured in that they reach a few microGauss in the 
most prominent structures such as galaxy clusters, but 
is < 1CT 11 G in the voids. We note that the EGMF in 
this LSS simulation arc relatively extended compared to 
other simulations [4lj and current observations do not 
allow to distinguish between such different EGMF sce- 
narios. The EGMF obtained in Ref. |4l| is closer to the 
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FIG. 3: Cross-section through the LSS simulation box from 
Refs. [1^, [4(| used in this study. The observer position is 
marked as a black cross near the center. The approximate 
positions of the ten sources in the discrete source realization 
investigated below are shown as black diamonds. Note that 
the simulation box including EGMF and sources is periodi- 
cally repeated in all three directions. Upper panel: Baryon 
density. Lower panel: magnetic field strength for the EGMF 
considered. The field polarizations are not shown and are in 
general a combination of smooth and stochastic components. 



case of negligible EGMF that we also study here. 

Since the LSS simulation covers a volume of only 
~ (75Mpc) 3 , the structures and fields in the simulation 
box are periodically repeated in all three directions. In 
this way, no evolution of the structure and fields with red- 
shift is taken into account. This is sufficient within the 
uncertainties in such scenarios since at distances larger 
than the GZK distance, spectra and angular distributions 
are insensitive to the detailed structure. The observer is 
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chosen in a low magnetic field region of the simulation 
box, resembling Earths actual environment, as can be 
seen in Fig. [3l 

io~ 5 1 — 1 — 1 — 1 — i 1 — i 
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FIG. 4: Upper panel: Fluxes of primary nucleons (black) 
and secondary photons (red) for continuous sources following 
the baryon density of the LSS simulation in Refs. 39, 40] and 
injecting an E~ 2 ' 6 proton spectrum up to 10 21 eV and redshift 
z = 3 with no source evolution, m — 0. Solid line assume the 
EGMF from these LSS simulations, whereas dashed lines are 
for negligible EGMF. Data as in Fig. [1] Lower panel: Photon 
to cosmic ray flux ratio for this case. 

Fig. H] shows results for structured, but continuous 
sources whose density is assumed proportional to the 
baryon density in the LSS simulation of which a cross- 
section is shown in Fig. [3l This time, the case with 
negligible EGMF is compared with the EGMF obtained 
from the LSS simulations, as in Refs. (38|. The spectra 
are hardly different from the homogeneous case shown in 
Fig. [TJ This scenario therefore corresponds to the uni- 
versal spectrum discussed in Ref . . The proton flux is 
only slightly suppressed by the EGMF at a few 10 18 eV. 
The neutral fraction, shown in the lower panel of Fig. [4] is 
also quite insensitive to both the continuous source dis- 
tribution and the EGMF. It is at the lower end of the 



range obtained in Ref. [2iJ [3(| . 

In the following we define the photon to cosmic ray flux 
ratio as relative to the total observed cosmic ray flux be- 
low 3 x 10 19 eV and relative to the simulated nucleon flux 
above 3 x 10 19 eV. This photon to cosmic ray flux ratio, 
also shown in Fig. |4j in fact broadly follows the attenua- 
tion length of the EM cascade [3l| . This is to be expected 
as long as the primary flux does not exhibit a strong 
break between the energy E at which the secondary to 
primary ratio is considered and the energies E' > E of 
the primaries mostly responsible for the secondary flux 
at energy E. Below we will see that below ~ 3 x 10 18 eV 
and above ~ 10 20 eV, the photon to nucleon ratio indeed 
depends little on source distributions and magnetization 
and mostly follows the behavior of the EM cascade at- 
tenuation length. Since in the combined CMB and URB 
this length scale increases roughly linearly with energy 
between ~ 10 20 eV and ~ 10 21 eV jsU, this leads to a 
rise in the photon to nucleon ratio up to the order of 1 
percent in this energy range. Unfortunately, due to the 
sharp drop of the primary cosmic ray flux expected at 
these energies, a photon to nucleon ratio of order percent 
may not be easier to detect than an order 10~ 4 ratio at 
lower energies. 

We now consider a specific realization with discrete 
sources of density n s ~ 3 x 10 _5 Mpc" 3 . As Fig.[3]shows, 
the sources in this realization are concentrated around 
the most prominent galaxy cluster. The spectra result- 
ing in this scenario are shown in Fig. upper panel. 
Both above the GZK energy around 10 20 eV and below 
~ 3 x 10 18 EeV the primary spectra are not universal 
any longer. This is understandable because the aver- 
age source distance ~ 30 Mpc is now comparable to the 
GZK distance. Further, a cosmic ray of energy E in a re- 
gion permeated by magnetic fields B has a Larmor radius 
rL ~ 1 (E /EeV)(B / fiG)" 1 kpc and during the lifetime of 
the Universe, t u ~ 10 Gyr, diffuses a distance 

/ p \ !/ 2 / r> \ -1/2 
W~(^ u ) 1/2 ~2(— ) (-) Mpc, 

(1) 

as long as this length scale is smaller than the spatial 
extent of the field. Fig. [3] shows that the structured 
EGMF reaches values of B ~ 0.1 /zG over several Mpc 
around the sources in the LSS simulation used here. Be- 
low ~ 3 x 10 18 eV, the diffusion length then also becomes 
smaller than the average source distance 

4^3.46 Mpc. (2) 

This leads to a suppression of the nucleon flux due to 
a partial confinement of UHECRs which was discussed 
qualitatively in Ref. [ll[. Note that this is not quite the 
same as the magnetic horizon effect discussed in Refs @, 
HH E E3] where typically uniform EGMFs were assumed 
which is less realistic in a structured universe. 

Above the threshold for pair production by photons 
with the CMB at ~ 10 15 eV and up to ~ 10 20 eV, the 
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FIG. 5: Same as Fig.[4]for the discrete source realization with 
10 sources in the LSS simulation box, see Fig. [3j injecting an 
E~ 2S proton spectrum up to 10 21 eV and redshift z = 3 with 
no source evolution, m = 0. This corresponds to a source 
density of ~ 3 x 10~ 5 Mpc -3 . A gain, solid line assume the 
EGMF from the LSS simulations, whereas dashed lines are 
for negligible EGMF, and no source evolution is assumed. 



photon attenuation length is smaller than 10 Mpc [3l| . 
Photons in this energy range thus have to be produced 
within < 10 Mpc from the observer. Above ~ 3 x 10 18 eV 
7— rays are mostly produced by nucleons above the GZK 
threshold whose attenuation length is smaller than the 
typical source distance Eq. This explains why the 
photon fraction of the UHECR flux between ~ 3 x 10 18 eV 
and ~ 10 20 eV decreases with decreasing discrete source 
density. Furthermore, magnetic fields > 10 -10 G effec- 
tively block cascade development and more EM energy is 
channeled into synchrotron radiation, ending up at GeV- 
TeV energies. This explains why the photon to charged 
cosmic ray ratio also decreases somewhat with increasing 
EGMF. Both tendencies are seen by comparing Figs. 2] 
and[5j lower panels. 

In contrast, below ~ 3 x 10 18 eV, a good fraction of the 
photons are produced by pair production of protons at 
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FIG. 6: Same as Fig. [5] for the discrete source realization, 
but now the spectra resulting for negligible EGMF (dashed 
lines) and with the EGMF (solid lines) are compared for dif- 
ferent injection parameters, namely m = 0, E~ 2 ' 7 for the case 
without EGMF and m = 3, £T 2 ' 7 for the case with EGMF. 



energies below the GZK threshold. These protons have 
attenuation lengths larger than the typical source dis- 
tance Eq. ([2]). Furthermore, at these energies the EM 
attenuation length decreases to below a few Mpc and the 
cascade flux thus depends mostly on the environment 
of the observer within a few Mpc where the EGMF is 
negligible, see Fig. [3] These two facts explain why be- 
low ~ 3 x 10 18 eV the photon to cosmic ray flux ratio 
is relatively insensitive to discrete source densities and 
magnetic environments, as confirmed by Figs. [4] and [5j 
lower panels. In contrast, this ratio does decrease with 
increasing magnetic fields around the observer if these 
fields are > 10 -10 G, as we will see below for an EGMF 
with homogeneous statistical properties. 

In Figs.[5]and[7]we compare the spectra in the scenarios 
with and without EGMF for different source evolutions, 
namely m = 3 and m — 0, respectively, for an injec- 
tion spectrum oc E~ 2 ' 7 . In the case with EGMF, even 
for strong redshift evolution, a new low-energy, presum- 



6 




FIG. 7: Same as Fig. [|J] upper panel, but with the energies 
measured by HiRes and AGASA multiplied by a factor of 1.2 
and 0.9, respectively, and showing the differential nucleon flux 
multiplied by E s . 



ably galactic component becomes necessary below a few 
10 18 eV. Furthermore, in this case the neutral fractions 
shown in the lower panel of Fig. [6] extend considerably 
below the lower limit obtained in Ref. 0, [3(| which did 
not consider structured sources and EGMF. 

Finally, in order to elucidate the role of the struc- 
ture of the EGMF, we consider an EGMF with homoge- 
neous properties: Its Fourier transform follows a power 
law spectrum, |B(k)| 2 oc fc~ 2 5 for (75 Mpc) -1 < k < 
(0.6 Mpc) -1 , with an r.m.s. strength of 10~ 9 G. This 
mimicks the scenarios considered in Ref. [ll[ of an EGMF 
with a coherence length ~ 1 Mpc. The field strength is at 
the high end of realistic values, but still consistent with 
current upper limits [HI]. Fig. [5] compares the resulting 
spectra with the case of negligble EGMF. Comparing the 
upper panels of Figs. Eland Fig. [51 we see that even in this 
relatively strong EGMF with homogeneous properties, 
the proton flux suppression due to the magnetic horizon 
effect is considerably more modest than the confinement 
effect for the structured EGMF envisaged here. 

In the lower panel of Fig. [5] we see that an EGMF of 
strength > 10 _10 G with homogeneous properties tends 
to suppress the photon to cosmic ray flux ratio also at 
energies below — 3 x 10 18 cV. as was found already in 
Ref. [3] and also seen in Ref. [lij [3(3]. As remarked 
above, this is to be expected because the EM cascade de- 
velopment at these energies mostly depends on the mag- 
netization of the environment of the observer within a 
few Mpc. 

Note that the injection spectra needed in the scenar- 
ios discussed here where the cross over to extra-galactic 
cosmic rays occurs at a few 10 17 eV is much steeper than 
predicted by most acceleration scenarios, but could be 
explained as an effective spectrum obtained by averaging 
over sources with different maximal energies and harder 
individual spectra 45 1. 
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FIG. 8: As in Figs. [5] and [6] for the discrete source realization 
injecting an E~ 2 ' 6 proton spectrum in the absence of redshift 
evolution, m — 0, but now the spectra resulting for negligible 
EGMF (dashed lines) are compared with the spectra for a 
stochastic EGMF of strength 10 -9 G and effective coherence 
length ~ 1 Mpc (solid line), see text for details. 



The proton injection power required for the scenar- 
ios shown in Figs. [5][7] is ~ 7 x 10 37 ergs~ 1 Mpc~ 3 for 
negligible EGMF and about a factor 10 higher, ~ 6 x 
10 38 erg s -1 Mpc -3 , in the presence of the EGMF consid- 
ered here. The number for negligible fields is consistent 
with the power obtained in Ref. [421 ] which only consid- 
ered relatively weak, unstructured EGMF. The orders 
of magnitude are consistent with acceleration in active 
galactic nuclei, for example. 

All the scenarios discussed here are consistent with 
present data in terms of large scale isotropy and auto- 
correlations. However, an increase of the current world 
exposure by a factor of a few predicts significant large 
scale anisotropy as well as small-scale clustering, at least 
in the case of modest EGMF. Absence of such effects 
would hint to a relatively strong EGMF, as demonstrated 
in Fig. [SI or to a considerably higher source density. 

We also found that the diffuse cosmogenic neutrino 
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FIG. 9: Autocorrelation function of UHECR in the discrete 
source scenarios shown in Figs. [5] and [6] assuming 10 3 events 
observed above 40 EeV, a factor 5 — 10 above the current 
world data set. An angular resolution of 1° was assumed. 
Upper panel: No EGMF. Lower panel: Assuming the EGMF 
from the LSS simulation. 



flux depends mostly on the source evolution, but little 
on the source distribution or EGMF. This is understand- 
able because the neutrino energy attenuation length is 
basically the Hubble radius and thus always larger than 
all other scales involved. Sources up to cosmological dis- 
tances thus contribute to the diffuse neutrino flux. Since 
neutrinos are produced by pionproduction of nucleons 
off the CMB and IR backgrounds, for a given UHECR 
injection spectrum, the fraction of the UHECR energy 
transformed to neutrinos should, therefore, be rather in- 
dependent of source distribution or EGMF. Our neutrino 
fluxes are comparable to other calculations in the liter- 
ature, for example Ref. AG]. A possible enhancement of 
cosmogenic neutrino fluxes from the source environment 
could come from an enhancement of the IR background 
around sources such as galaxy clusters, an effect we have 
not taken into account in the present work. One case in 
which this possibility was considered [JtJ does, however, 



not give much higher diffuse neutrino fluxes when the 
constraint from the primary UHECR flux is taken into 
account. 



IV. UNCERTAINTIES IN THE PHOTON TO 
CHARGED COSMIC RAY FLUX RATIO 

If the actual but poorly known URB is larger than the 
minimal URB we used in the simulations presented in 
the previous section, photon absorption will be stronger, 
further reducing the photon to cosmic ray fraction. 

A significant fraction of heavier nuclei of atomic num- 
ber A in the primary UHECR flux is likely to reduce the 
photon to charged cosmic ray flux ratio above ~ 10 19 eV: 
Pion-production is shifted down by a factor ~ A in en- 
ergy, falling below ~ 10 19 eV already for moderate A. An 
additional source of 7— rays is from photo-spallation of 
nuclei |48[ . But these photons have an energy of ~ Ta x 
a few MeV which for Lorentz factors < 10 11 is below 
~ 10 18 eV. We recall in this context that the cosmogenic 
neutrino flux is also believed to considerably depend on 
primary cosmic ray composition (49l . [Hoi . l5lT ] . 

On the other hand, additional photons could be pro- 
duced in nucleon interactions with an IR background that 
could be enhanced around sources such as galaxy clus- 
ters, or by nucleon interactions with the ambient baryon 
gas in the sources. Both these possibilities have not been 
taken into account in the present work. However, we do 
not expect that these effects could significantly increase 
the photon fraction above ~ 10 15 eV for the following 
reasons: First, as discussed in Sect. Ill above, since in 
this energy range, the EM cascade attenuation length 
is smaller than or comparable to the typical source dis- 
tance, the photon fraction depends mostly on EM sec- 
ondary production outside the sources. Second, the pri- 
mary cosmic ray flux is not expected to be significantly 
modified by these processes, since only about 1% of its 
energy is expected to be converted by interactions with 
an enhanced IR background within galaxy clusters [47| . 
and since the optical depth for pp interactions is less than 
unity in galaxy clusters [52| . 



V. CONCLUSIONS 

We have performed simulations following ultra-high 
energy nucleon trajectories above 10 17 eV and compared 
scenarios with continuous and discrete source distribu- 
tions as well as the case of negligible and highly struc- 
tured extragalactic magnetic fields reaching microGauss 
levels in galaxy clusters. Powerful ultra-high energy cos- 
mic ray sources likely have densities of order 10~ 5 Mpc -3 . 
We find two considerably non-universal effects that de- 
pend on source distributions and magnetic fields: We 
found that the primary extra-galactic cosmic ray flux 
can become strongly suppressed below a few 10 18 eV if 
the sources are immersed in highly structured magnetic 
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fields. This is due to the fact that cosmic ray primaries 
start to be magnetically confined and do not reach the ob- 
server any more at these energies. This effect can indeed 
be considerably stronger than for stochastic magnetic 
fields with homogeneous properties and r.m.s. strength 
< 10~ 9 G. We also found that the secondary photon to 
primary cosmic ray flux ratio between ~ 3 x 10 18 eV and 
~ 10 20 eV is of the order 10~ 4 and decreases with de- 
creasing source density and increasing magnetization. In 
principle at least, this ratio could therefore serve as an 
independent measure of these poorly known quantities. 
In contrast, the photon to cosmic ray flux ratio is rather 
insensitive to source distributions and magnetic fields at 
lower energies as well as above ~ 10 20 eV where it can 
reach the percent level. We also pointed out that a sig- 
nificant contribution of nuclei heavier than hydrogen to 
the primary cosmic ray flux tends to reduce the 7— ray to 
charged cosmic ray ratio further. Also, additional pho- 
ton production within the sources is unlikely to increase 



this ratio significantly. As a consequence, in acceleration 
scenarios for the origin of highest energy cosmic rays the 
fraction of secondary photons may be difficult to detect 
even for experiments such as Pierre Auger. In contrast, 
the cosmogenic neutrino flux does not significantly de- 
pend on source density and magnetization, but mostly 
on the redshift evolution of the sources. 
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